Magnetic vortex echoes: application to the study of arrays of magnetic nanostructures 
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We propose the use of the gyrotropic motion of vortex cores in nanomagnets to produce a magnetic 
echo, analogous to the spin echo in NMR. This echo occurs when an array of nanomagnets, e.g., 
nanodisks, is magnetized with an in-plane (xy) field, and after a time r a field pulse inverts the 
core magnetization; the echo is a peak in M xy at t — It. Its relaxation times depend on the 
inhomogeneity, on the interaction between the nanodots and on the Gilbert damping constant a. 
Its feasibility is demonstrated using micromagnetic simulation. To illustrate an application of the 
echoes, we have determined the inhomogeneity and measured the magnetic interaction in an array 
of nanodisks separated by a distance d, finding a d~ n dependence, with n « 4. 



PACS numbers: 75.70.Kw,75.78.Cd,62.23.Eg,76.60.Lz 

The interest in magnetic vortices and their properties 
and applications has grown steadily in the last years [lj- 
4]. Vortices have been observed, for example, in disks 
and ellipses having sub-micron dimensions [5j. More re- 
cently, the question of the intensity of the coupling be- 
tween neighbor disks with magnetic vortex structures has 
attracted an increasing interest 

At the vortex core the magnetization points perpendic- 
ularly to the plane: this characterizes its polarity, p = +1 
for the +z direction and p — —1 for — z. The direc- 
tion of the moments in the vortex defines the circulation: 
c = — 1 for clockwise (CW) direction, and +1 for CCW. 
If removed from the equilibrium position at the center 
of the nanodisk by, for example, an in-plane field, and 
then left to relax, a vortex core will perform a gyrotropic 
motion, with angular frequency to, given for thin dis ks [3| 
by u G w (20/9)jM s /3 (J3 = h/R is the aspect ratio) \v$. 

We propose in this paper that, manipulating the dy- 
namic properties of the vortex in an analogous way as it 
is done in Nuclear Magnetic Resonance (NMR), a new 
phenomenon results, the magnetic vortex echo (MVE), 
similar to the spin echo observed in NMRfllJ. This new 
echo may provide information on fundamental properties 
of arrays of nanodisks, e.g., their inhomogeneity and in- 
teractions. Despite the fact that applications of vortices 
necessarily involve arrays, most of the recent publications 
deal with the analysis of individual nanodisks or arrays 
with a few elements. Therefore, the possibility of char- 
acterizing large arrays is of much interest. In this paper 
we have examined the motion of vortex cores in an ar- 
ray of nanometric disks under the influence of a pulsed 
magnetic field, using micromagnetic simulation. 

Let us consider an array of nanodisks where the vor- 
tex cores precess with a distribution of angular frequen- 
cies centered on ojq, of width Aw, arising from any type 
of inhomogeneity (see note [10]). We assume that the 
frequencies vary continuously, and have a Gaussian dis- 
tribution P(u>) with mean square deviation Aco. To sim- 
plify we can assume that the polarization of every vortex 



FIG. 1. (Color online) Diagram showing the formation of 
magnetic vortex echoes; the disks are described from a refer- 
ence frame that turns with the average translation frequency 
cjo: a) disks with in-plane magnetizations M,: along the same 
direction (defined by the white arrows); b) after a time r 
the disks on the left, center and right have turned with fre- 
quencies, respectively, lower, higher and equal to ojo; c) the 
polarities of the vortex cores are reversed, and the uii of the 
vortex cores (and of the Mi) change sign, and d) after a sec- 
ond interval r the cores (and Mi) are again aligned, creating 
the echo. 



is the same: pi = +1. This is not necessary for our argu- 
ment, but, if required, the system can be prepared; see 
ref. [12| and the references therein. 

Since the direction of rotation of the magnetic vortex 
cores after removal of the in-plane field is defined exclu- 
sively by p, all the cores will turn in the same direction; 
as the vortex core turns, the in-plane magnetization of 
the nanodot also turns. 
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one can derive the total in-plane magnetization: 



10 20 30 40 50 60 70 80 90 
t(ns) 

FIG. 2. (Color online) Micromagnetic simulation of magnetic 
vortex echoes, for 100 nanodisks, with d = infinity, and a) 
a = 10 nm, r = 30 ns, a = 0; b) a — 20 nm, r = 10 ns 
and r — 40ns (two pulses), and a = 0.001; c) a = 20 nm, 
r = 20 ns, a = 0; d) a = 10 nm, r = 20 ns, a = 0.005. The 
inversion pulses (B z = — 300 mT) are also shown (in red). 
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FIG. 3. (Color online) Variation of the inverse of the re- 
laxation times T2 (diamonds) obtained by fitting the curves 
of echo intensity versus time interval r to Moexp (— T/T2), 
as a function of a, for simulations made with D — 250 nm, 
a — 10 nm, separation infinite; the continuous line is a linear 
least squares fit. 



Let us assume that all the vortex cores have been dis- 
placed from their equilibrium positions along the positive 
x axis[l3j. by a field B. The total in- plane magnetiza- 
tion (that is perpendicular to the displacement of the 
core) will point along the y axis, therefore forming an 
angle 9 = at t = 0. Using the approach employed in 
the description of magnetic resonance (e.g., see [lj, [lj|), 



My(t) = M v (0) 
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■ cos(ujt)duj (1) 



an integral that is[16j| the Fourier transform of the func- 
tion P(u); using T* = 1/Aw: 

M y (t) = M y (0) exp ("^) ( 2 ) 

This result shows that the total magnetization tends to 
zero, as the different contributions to M y (t) get gradu- 
ally out of phase. This decay is analogous to the free 
induction decay (FID) in NMR; its characteristic time is 
T 2 * = 1/Aw. 

After a time r, the angle rotated by each vortex core 
will be wt; if at t — r we invert the polarities of the vor- 
tices in the array, using an appropriate pulse, the motion 
of the cores will change direction (i.e., to — > —lj), and one 
obtains: 



My(t-T)=My(0) 



, AwV27r 

The magnetization at a time t > r is then: 



cos[w(r — t)]du> (3) 



M y (t) = M v (0) e 
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e ( T 2 ] cos(uj t) (4) 



This means that the magnetization component M y (t) in- 
creases for t < t < 2t, reaching a maximum at a time 
t = 2t: this maximum is the magnetic vortex echo, anal- 
ogous to the spin echo observed in magnetic resonance, 
which has important applications in NMR, including in 
Magnetic Resonance Imaging (MRI) fLU Il5| (Fig. [J). In 
the case of the NMR spin echo, the maximum arises from 
the refocusing of the in-plane components of the nuclear 
magnetization. 

In Eq. 2]we have included a relaxation term contain- 
ing the time constant T2 - also occurring in NMR -, to 
account for a possible decay of the echo amplitude with 
time; its justification will be given below. 

In the array of nanodisks, there will be in principle two 
contributions to the defocusing of the magnetization, i.e., 
two mechanisms for the loss of in-plane magnetization 
memory: 1) the spread in values of /3 and H (see note 
[10]), producing an angular frequency broadening term 
Au, and 2) irreversible processes that are characterized 
by a relaxation time T 2 : thus 1/T| = Acu + 1/T 2 . 

The second contribution is the homogeneous term 
whose inverse, T%, is the magnetic vortex transverse re- 
laxation time, analogous to the spin transverse relaxation 
time (or spin-spin relaxation time) T 2 in magnetic reso- 
nance. The irreversible processes include a) the interac- 
tion between the disks, which amounts to random mag- 
netic fields that will increase or decrease w of a given 
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disk, producing a frequency spread of width Aw' = l/T^, 
and b) the loss in magnetization (of rate arising 
from the energy dissipation related to the Gilbert damp- 
ing constant a. Identifying T n to the NMR longitudinal 
relaxation time T x , one has Q: 1/T 2 = l/T£ + 1/2T Q . 
Therefore the relaxation rate 1/T 2 * is given by: 

— = Aw + — = Aw + — + — (5) 

J 2 ±2 1 2 LI a 

1 /T 2 * is therefore the total relaxation rate of the in-plane 
magnetization, composed of a) Aw, the inhomogeneity 
term, and b) I/T2, the sum of all the other contributions, 
containing 1/T 2 , due to the interaction between the disks, 
and 1/T a , the rate of energy decay. The vortex cores 
will reach the equilibrium position at r — after a time 
t ~ T a , therefore there will be no echo for 2r 3> T a . 

The vortex echo maximum at t = 2r, from Eq. 2J is 
My(2r) oc exp(~ T/T2); one should therefore note that 
the maximum magnetization recovered at a time 2r de- 
creases exponentially with Ta, i.e., this maximum is only 
affected by the homogeneous part of the total decay rate 
given by Eq. [5] In other words, the vortex echo cancels 
the loss in M y due to the inhomogeneity Aw, but it does 
not cancel the decrease in M y due to the interaction be- 
tween the nanodisks (the homogeneous relaxation term 
l/Tg), or due to the energy dissipation (term l/2T a ). 

Note also that if one attempted to estimate the inho- 
mogeneity of an array of nanodots using another method, 
for example, measuring the linewidth of a FMR spec- 
trum, one would have the contribution of this inhomo- 
geneity together with the other terms that appear in Eq. 
arising from interaction between the dots and from 
the damping. On the other hand, measuring the vor- 
tex echo it would be possible to separate the intrinsic 
inhomogeneity from these contributions, since Ta can be 
measured separately, independently of the term Aw. T2 
can be measured by determining the decay of the echo 
amplitude for different values of the interval r. 

The Fourier transform of either the vortex free induc- 
tion decay or the time dependence of the echo M y (t) gives 
the distribution of gyrotropic frequencies P{u>). 

For the experimental study of vortex echoes, the se- 
quence of preparation (at t=0) and inversion fields (at 
t=r) should of course be repeated periodically, with a 
period T ^> T a . As in pulse NMR, this will produce 
echoes on every cycle, improving the S/N ratio of the 
measured signals. Also note that the time T 2 * can be ob- 
tained either from the initial decay (FID, Eq. [2|) or from 
the echo (Eq. |4j, but T2 can only be obtained from the 
MVE. 

In order to demonstrate the MVE, we have performed 
micromagnetic simulations of an assembly of 100 mag- 
netic nanodisks employing the OOMMF codejl?]. The 
simulated system was a square array of 10 x 10 disks, 
thickness 20 nm, with distance d from center to center. 
In order to simulate the inhomogeneity of the system, 



we have introduced a Gaussian distribution of diame- 
ters, centered on 250 nm and mean square deviation a; 
a = 10 nm corresponds to Aw w 1.6 x 10 8 s _1 . The 
disks were placed at random on the square lattice. The 
initial state of the disks (p = +1 and c = — 1) was pre- 
pared by applying an in-plane field of 25 mT; the po- 
larity was inverted with a Gaussian pulse of amplitude 
B z = -300 mT, with width 100 ps. The results for the 
case d = 00 were simulations made on the disks one at a 
time, adding the individual magnetic moments fii(t). 

We have successfully demonstrated the occurrence of 
the magnetic vortex phenomenon, and have shown its 
potential as a characterization technique. The simula- 
tions have confirmed the occurrence of the echoes at the 
expected times (t = 2r). For different values of a, the 
T 2 * time, and consequently the duration of the FID and 
the width of the echo are modified (Fig. [2^, [2};) ; increas- 
ing a results in a faster decay of the echo intensity as a 
function of time (Fig. [2^, EJi). We have also obtained 
multiple echoes, by exciting the system with two pulses 

(Fig. Eb)Q. 

Fig. |3] shows the dependence of T2 on a for a = 10 nm; 
essentially the same result is obtained for a — 20 nm, 
since T 2 does not depend on Aw (Eq. [5]). Taking a linear 
approximation, 1/T 2 = Aa, and since for d ^infinity 
there is no interaction between the disks, I/T2 = 1/2T Q , 
and therefore: 

-L = 2Aa (6) 

From the least squares fit (Fig. A = 1.6 x 10 10 s _1 . 
This relation can be used to determine experimentally a, 
measuring T 2 with vortex echoes, for an array of well- 
separated disks. 

Recently some workers have analyzed the important 
problem of the interaction between disks exhibiting mag- 
netic vortices, obtaining that it varies with a d~ n depen- 
dence: Vogel and co-workers [6j, using FMR, obtained 
for a 4 x 300 array a dependence of the form d~ 6 , the 
same found by Sugimoto et al. 0] using a pair of disks 
excited with rf current. Jung et al. studying a pair of 
nanodisks with time-resolved X-ray spectroscopy, found 
n = 3.91 ± 0.07 and Sukhostavets et al. Q, also for 
a pair of disks, in this case studied by micromagnetic 
simulation, obtained n = 3.2 and 3.7 for the x and y 
interaction terms, respectively. 

As a first approximation one can derive the depen- 
dence of the contribution to 1 /T 2 * related to the distance 
between the disks as: 

T* =B + CdT n (7) 

From our simulations, and using Eq. [7] we found, from 
the best fit, that this interaction varies as d~ n , with n = 
3.9 ± 0.1, in a good agreement with Q and reasonable 
agreement with Sukhostavets et al.Q. 
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FIG. 4. (Color online) Variation of the relaxation times 
TJ versus d~ 4 for an array of 10 x 10 magnetic nanodisks 
with a distribution of diameters centered on D — 250 nm 
(<r = 10 nm), damping constant a = 0.001 and separation 
d between their centers; the continuous line is a linear least 
squares fit. The inset shows a vortex echo simulation for the 
array, with d = 500 nm, r = 30 ns, a = 0.001. 



Determining XjT^ has allowed us to obtain the inten- 
sity of the interaction between the disks as a function 
of separation d between them. Substituting Eq. [5] and 
Eq. [7] in Eq. O we can obtain the expression for the 
interaction as a function of d: 



1 

1 2 



1 



B + Cd- 



Aa - Auj ps — - Aa - Acj; (8) 
1^1 



(approximation valid for d small). In Fig. [4] we show the 
results of the simulations with a = 10 nm and a = 0.001. 
Assuming n = 4 and making a linear squares fit, we 
obtained B = 6.15 x 10~ 9 s, C = -4.03 x 10~ 35 s m 4 . 

A new phenomenon, the magnetic vortex echo, anal- 
ogous to the NMR spin echo, is proposed and demon- 
strated here through micromagnetic simulation. Appli- 
cations of the magnetic vortex echo includes the measure- 
ment of the inhomogeneity, such as, distribution of di- 
mensions, aspect ratios, defects, and perpendicular mag- 
netic fields and so on, in a planar array of nanodisks or 
ellipses; it may be used to study arrays of nanowires or 
nanopillars containing thin layers of magnetic material. 
These properties cannot be obtained directly, for exam- 
ple, from the linewidth of FMR absorption. 

The MVE is a tool that can be used to evaluate the 
interaction between the elements of a large array of nano- 
magnets with vortex ground states. It can also be used 
to determine the Gilbert damping constant a in these 



systems. 
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